Hadron multiplicity in pp and AA collisions at LHC from the Color Glass Condensate 
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We provide quantitative predictions for the rapidity, centrality and energy dependencies of inclu- 
sive charged-hadron productions for the forthcoming LHC measurements in nucleus-nucleus colli- 
sions based on the idea of gluon saturation in the color-glass condensate framework. Our formulation 
gives very good descriptions of the first data from the LHC for the inclusive charged-hadron pro- 
duction in proton-proton collisions, the deep inelastic scattering at HERA at small Bjorken-a;, and 
the hadron multiplicities in nucleus-nucleus collisions at RHIC. 



I. INTRODUCTION 

It has been vifidely discussed that QCD predicts at 
high-energy formation of a nevi' state of matter, the so- 
called color-glass condensate (CGC) [l|. In this state, 
the density of quarks and gluons p with transverse mo- 
menta less that Qs reaches a high value, namely p ex 
l/as{Qs) S> 1 where as is the strong coupling constant 
and Qs is a new momentum scale (saturation momen- 
tum) that increases with energy. Therefore, as (Qs) ^ 
1 and this fact allows us to treat this system on solid the- 
oretical basis. One of the most characteristic and qual- 
itative consequences of the CGC is the emergence of a 
new mechanism for hadron production at high energy. 
In this approach, the process of secondary hadron pro- 
duction goes in two stages: production of gluon mini-jet 
with a typical transverse momentum Qs; and the decay 
of gluon mini-jets into hadrons. 

The CGC picture for the hadron production has passed 
two critical tests. First, it explains the main features of 
hadron multiplicity in heavy ion-ion collisions at RHIC 
(KLN papers 2]); and it predicts the inclusive hadron 
production in proton-proton (pp) collisions 'Sj at the LHC 
which was recently confirmed experimentally by the CMS 
collaboration at ^/s = 7 TeV [j] and the lower LHC en- 
ergies p-Hl' see also Ref. Q- The ion-ion collisions at 
the LHC will lead to a crucial test of the CGC approach. 
Indeed, the KLN success ^ in the description of RHIC 
data shows that we are dealing with scatterings of dense 
partonic systems even at rather low energies. At higher 
energies we expect that the density of scattering systems 
increases resulting in more transparent and unambiguous 
evidence for the creation of the CGC. 

In this paper we wish to extend our approach [3] to 
ion-ion (A A) collisions at the LHC. As it was discussed 
in Ref. Q, we improve the KLN approach in various 
ways incorporating the unintegrated gluon density that 
describes HERA data at small Bjorken-x. We also in- 
clude our knowledge on proton-proton scatterings at the 
LHC to make our prediction more reliable. In particular, 
we re-calculate the saturation momentum for nuclei from 
the saturation momentum in the proton target obtained 
from DIS data at HERA and confront it with the RHIC's 



gold-gold data. On the contrary, in the KLN approach 
[9j , the LHC saturation momentum was found via an ex- 
trapolation of the energy dependence of the saturation 
scale at RHIC in the BFKL region. 

In the next section, we discuss the /cT-factorization and 
our main formulation for the inclusive hadron produc- 
tion in pp and AA collisions. In particular, we introduce 
the impact-parameter dependent saturation scale for a 
proton and a nuclear target. Section HI is devoted to 
comparison with the experimental data and to discus- 
sion of various predictions for the LHC energies. Finally, 
we conclude in Sec. IV. 



II. MAIN FORMULATION 

The gluon jet production in nucleus-nucleus collisions 
can be described by /cT-factorization given by flOll , 
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where a;i,2 = {pt / V^)e^^ , pt and y are the transverse- 
momentum and rapidity of the produced gluon mini-jet. 
(f>^{xi; kr) denotes the unintegrated gluon density and is 
the probability to find a gluon that carries Xi fraction of 
energy with kx transverse momentum in the projectile 
(or target) nucleus A. We defined Cp = {N^ - l)/2A^c 
where Nc denotes the number of colors. 

The /cT-factorization has been proven [10| for the scat- 
tering of a diluted system of parton with a dense one. 
Such a process can be characterized by two hard scales: 
the transverse momentum of produced particle px and 
a saturation scale which are both larger than the soft 
interaction scale p. The ion-ion scattering is a typical 
example in which we have three scales: pt and two sat- 
uration scales for the projectile and the target. The fc^- 
factorization might be then violated only at low xi and 
X2, and when pt is smaller than both saturation scales. 
Outside of this kinematic region we are actually dealing 
with scatterings of diluted-dense system since one of the 
saturation scales is small. For the case of scatterings of 
dense-dense system of partons, the /cT-factorization has 
not yet been proven. Nevertheless, the success of the 



KLN approach 2] which is based on the /cT-factorization, 
in description of the experimental data at RHIC for gold- 
gold collisions suggests that the fcT-factorization is cur- 
rently the best tool that we have at our disposal for the 
processes of ion-ion scatterings. 

The relation between the unintegrated gluon density 
(/)^ and the color dipole-nucleus forward scattering am- 
plitude has been obtained in Ref. [13] • It reads as follows jy . (x- r-b) = < 



amplitude [l3|. This model describes both the HERA 
DIS data at small- a; 14] and the proton-proton LHC data 
^3] . Second, we replace the proton saturation momentum 
by that of the nucleus. The dipole-nucleus scattering am- 
plitude in our model is given by 
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NX {x,; rr; b) = 2Na (x^; tt; b) - N^ {x,-rT] b) , (3) 

where NA{xi;rT',b) is the dipole-nucleus forward scat- 
tering amplitude which satisfies the perturbative nonlin- 
ear small-x Balitsky-Kovchegov (BK) quantum evolution 
equation lll|. In the above, ry denotes the dipole trans- 
verse size and b is the impact parameter of the scattering. 
Substituting Eq. ([2]) into Eq. ([T]), and after analytically 
performing some integrals, we obtain |3l. llOl]. 
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where B is the impact parameter between the center of 
two nuclei, 6 and b^ = b — B are the impact parameter 
between the interacting nucleons with respect to the cen- 
ter of two nuclei. We will see below that we can neglect 
the impact parameter of the produced gluon-jet from the 
center of the nucleons. 

In the above, we extended the fcT-factorization given in 
Eq. ([T]) by introducing a running strong-coupling as ■ For 
the running as we employ the same prescription used in 
Ref. [3] for pp collisions. The saturation scale QAixi^2,b) 
depends on the a;i.2 and the impact parameter and will 
be introduced in the following. A given centrality bin cor- 
responds to ranges of the impact-parameter B G [Bi, B2] 
of the collisions. The procedure how to select the event 
with a fixed B and its relation to the centrality bins is 
well-known, see for example Refs. |9l[l^. 

Notice that the relation between the unintegrated 
gluon density and the forward dipole-nucleus amplitude 
Eqs. (I2I3P in the fcT-factorization Eq. ([T]) is not a simple 
Fourier transformation which is commonly used in liter- 
ature and also depends on the impact-parameter. The 
impact-parameter dependence in these equations is not 
trivial and should not be in principle assumed as an over- 
all factor. Using the general properties of high density 
QCD that the underlying physics depends only on the 
saturation scale we reconstruct the dipole-nucleus scat- 
tering amplitude in two steps. First, we choose the sat- 
uration model which effectively incorporates all known 
saturation properties [3| driven by the BK equation in- 
cluding the impact-parameter dependence of the dipole 
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where we defined Z — rQA{x;b), Y — ln(l/x) and 
K = x"{ls)/x'{ls) where x is the LO BFKL character- 
istic function. The parameters A and B are determined 
uniquely from the matching of Na and its logarithmic 
derivatives at Z = 2. The nucleus saturation scale is 
given by 



Ql{x;b) 
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In the above, the proton saturation scale Qp is defined 



as 



Qpix;b') = (^)' 



exp 
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and Ta {B) denotes the nuclear thickness. We use for the 
nuclear thickness the Wood-Saxon parametrization [l6|. 
Notice that for small Z < 2, the effective anomalous 
dimension jcr = 1 — 7s in the exponent in the upper line 
of Eq. ([5]) rises from the BK value towards the DGLAP 
value. The ansatz given in Eq. ([5]) for a proton target was 
first introduced in Ref. |15| . The parameters A, 7^, A^o, a^o 
and BcGC are obtained from a fit to the DIS data at low 
Biorken-x x < 0.01 with a very good x^/d.o.f. = 0.92 

The saturation scale Qa in nuclei is proportional to 
the density of partons in the transverse plane (see for 
example Refs. [ll, 0, 113]) and Eq. ([6]) takes this fact into 
account. In Eq. ([6]), b' is the impact parameter of the 
color dipole with respect to the center of the proton, b 
is the impact parameter of the nucleon with respect to 
the center of the nucleus and b — b' is, therefore, the 
position of the dipole with respect to the center of the 
nucleon. The nuclear saturation scale defined in Eq. (j6]) 
gives Q\ ~ QpA^^^, in agreement with the basic idea 
of saturation and the CGC picture [l|, Q, ll3- This is 
because Ta ~ A^'^ (where A is the effective mass num- 
ber of the nucleus in a given centrality) and the integral 
in Eq. ^ is approximately over the nucleon size. Note 
that since here the nucleon saturation scale depends on 
impact-parameter, it is not then correct to define the 
saturation scale on a nuclear target with a simple pref- 
actor A^/^ scaling. In Fig. [1] we show that the impact- 
parameter dependence of the saturation scale for proton 
and lead obtained from Eqs. (I6I7I) are quite different as we 
expected. Fig. [T] demonstrates the well-known features 
of nucleus scattering, namely the typical b' in Eq. ^ is 
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presence of a non-zero mini-jet mass, namely replacing 
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in a;i,X2 and also in the denominator 
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FIG. 1: The impact parameter dependence of the saturation 
scale for proton and lead at a; = 10~^ and x = 10~^. 



much smaller than the typical b ^ Ra {Ra is the nuclear 
radius) in the nuclear thickness. 

In order to simulate the behavior of gluon density at 
large a; — ?> 1, we product the unintegrated gluon density 
with (1 — x)"^ as prescribed by quark counting rules and 
the HERA data on DIS at large-x [1]. Notice that the 
contribution of (1 — x)* correction at ^/s = 5.5 TeV for 
AA collisions and rj < 4.5 is less than 4%, and for 4.5 < 
7^ < 6 is less than 13%. Therefore, the main contribution 
of the unintegrated gluon density at the LHC high-energy 
in the kinematic region considered here for both AA and 
also pp collisions [3J comes from the small- a: region where 
the saturation physics is important. For the importance 
of the saturation effects in pA collisions at the LHC, see 
Ref. fl^. 

The fcr-factorization Eqs. (|ll4p gives the cross-section 
of radiated gluon mini-jets with zero mass while what 
is actually measured experimentally is the distribution 
of final hadrons. We therefore should model the non- 
perturbative hadronization stage of gluon mini-jets. We 
employ the Local Parton-Hadron Duality principle [19| . 
namely the hadronization is a soft process and cannot 
change the direction of the emitted radiation. Hence 
the form of the rapidity distribution of the mini-jet and 
the produced hadron are different only with a numeri- 
cal factor C. It is well-known that the general assump- 
tion about hadronization leads to the appearance of mass 
of the mini-jet which is approximately on average equal 
to m|gj ~ 2fi{pT) [2i, l^l where fi is the scale of soft in- 
teraction. The mini-jet mass rrijet effectively incorpo- 
rates the nonperturbative soft pre-hadronization in the 
pseudo-rapidity space and can be approximately related 
to the saturation scale Q- Accordingly, one should also 
correct the kinematics every where in Eq. Q due to the 



of 1/px- Finally, in order to take account of the difference 
between rapidity y and the measured pseudo-rapidity ry, 
we employ the Jacobian transformation between y and 77 

i. 

Following Ref. Q, in the spirit of the geometrical- 
scaling property of the scattering amplitude, we obtain 
the charged-particle multiplicity distribution at a fixed 
centrality but various energies from the corresponding 
mini-jet cross-section Eq. ^ divided by the average area 

of interaction cr^ = M-nUr^^^^ — Mtt (b'^ + b'^_) . In a 

similar way, in order to obtain the charged-particle mul- 
tiplicity distribution at various centrality bins but a fixed 
energy, we divide the mini-jet cross-section integrated in 
ranges oi B G [Bi,B2] with the corresponding relative 
interaction area ag = A^7r(i?| — Bf). 



III. DISCUSSION AND PREDICTIONS 

We have only two unknown parameters in our model: 
the overall factor ^ and the soft-scale /x introduced 
in the definition of mini-jet mas. We also introduced 
a A'-factor which incorporates the discrepancy between 
the exact calculation with our formulation. These two 
phenomenological parameters are fixed at RHIC energy 
yi = 200 GeV for Au-Au - 6% centrality Then our 
results at other energies and various centralities can be 
considered as free-parameter predictions. The sensitivity 
of our results to various /i is shown in Fig. [5] (top) at 
y/s ~ 200 GeV. The preferred value of /j, in Pb-Pb col- 
lisions is approximately of order the current-quark mass 
and is smaller than the corresponding value of /x « m,r 
in pp collisions [3] . The decrease in the value of the soft- 
scale in a denser medium is in accordance with the notion 
of asymptotic deconfinement that the confinement radius 
increases with density. The preferred value of the soft- 
scale in AA collisions is very close to what was obtained 
in Ref. [201 ^r the value of Aqcd at moderate densities. 
Loosely-speaking, the smallness of p, for AA collisions 
compared to the pp collisions is due to the fact that the 
mini-jet mass reduces in a denser medium despite the 
fact that the saturation scale increases. We checked that 
for Pb-Pb and pp collisions a fixed mini-jet mass about 
0.14 — 0.2 and 0.4 — 0.5 GeV respectively, gives similar 
results to the case that the mini-jet mass runs with the 
saturation scale. As we already pointed out the mini- 
jet mass mimics the properties of the soft nonperturba- 
tive pre-hadronization stage and a different value for /i 
in AA collisions is simply due to the fact that the soft 
pre-hadronization stage in AA collisions is different from 
the corresponding one in pp collisions. 

The main sources of theoretical uncertainties in our 
approach stem from assuming a fixed /i (or a fixed 
rrijet) and the prefactor -^ for all energies, rapidi- 
ties and centralities. The experimental errors in the 
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FIG. 2: Top: The effect of the soft-scale /i is shown at ^/s — 
200 GeV. Lower: Pseudo-rapidity distribution of charged par- 
ticles produced in Au-Au and Pb-Pb central — 6% colli- 
sions at RHIC ys = 19.6, 130, 200 GeV and the LHC energies 
y/s = 2.75, 5.5 TeV. The band indicates less than 3% theoret- 
ical error coming from uncertainties related to normalization 
and modeling the mini-jet mass. The experimental data are 
from PHOBOS collaboration [H. 



data points taken for fixing these two unknown param- 
eters also induce some uncertainties. Moreover, we em- 
ployed a Glauber based approach to relate the centrality 
bins, number of participant and the cut in the impact- 
parameters [9|, uM- This formulation agreed with the 
RHIC data on number of participant within error bars, 
but we expect some uncertainties in the Glauber formu- 
lation mainly due to its over-simplicity. In the Glauber 
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FIG. 3: The pseudo-rapidity dependence at RHIC -^s = 200 
GeV (top) and the LHG y^ = 5.5 TeV (lower) at different 
centrality bins. The band indicates less than 3% theoretical 
errors. The experimental data are from PHOBOS collabora- 
tion fil. 



formalism one should also assume the value of inelastic 
nucleon-nucleon cross-section crj^'jf' (without diffractive 
components) from outset. Here we used the values ob- 
tained in Ref. 215: cr™'^' = 64.8,58.5,42,41,30 mb for 
^/i = 5500, 2750, 200, 130, 19.6 GeV, respectively. Notice 
that the nuclear saturation scale defined in Eq. ([6]) can 
be in principle different with exact one with extra factor. 
In order to estimate this uncertainty, we checked that a 
factor of 2 difference in the definition of Eq. ^ (namely 
Qa -^ "^Qa) will change our result less than 5% at the 
LHC. Overall, including all the above-mentioned possible 
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FIG. 4: Energy dependence of the charged hadrons multiplic- 
ity at midrapidity r; = in central collisions in pp and AA 
collisions. The theoretical curve Saturation model (LR) is our 
prediction. The band indicates less than 3% theoretical er- 
rors. The total theoretical uncertainties is less than 7%, see 
the text for the details. We also show the KLN prediction [0| 
with the same error band as ours. The experimental data are 
from [Hi,[23,[2l. 



theoretical uncertainties, we expect less than 7% theoret- 
ical error in our calculations at the LHC energies. The 
theoretical bar in Fig. l3l4l onlv show about 3% theoretical 
error. 

In Fig. [2] (lower), we show our predictions at lower 
RHIC energies ^s = 19.6 and 130 GeV in Au-Au colli- 
sions, and also for the LHC energies y/s — 2.75 and 5.5 
TeV in Pb-Pb collisions for — 6% centrality bin. 
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FIG. 5: The scaled pseudo-rapidity density as a function 
of number of participant Npar at midrapidity for Au-Au at 
^/s = 130, 200 GeV and for Pb-Pb at ^ = 5.5 TeV. The 
experimental data are from PHOBOS collaboration [2q |. 



In Fig. [3J we show the charged-particle multiplicity 
at various centrality bins for RHIC and the LHC high 
energy. As we already stressed we only used RHIC data 
at y/s = 200 GeV for — 6% centrality bin in order to 
fix two unknown parameter of our model. Therefore, the 
description of other RHIC data here can be considered 
as predictions. Over-all, our approach gives a very good 
description of RHIC multiplicity data as a function of 
rapidity, centrality, number of participant and energy, 
see Figs. UlSl 

In Fig. m we show the energy dependence of dNpp/drj, 
dNAA/dri and {2/Npar)dNAA/dr] at midrapidity 77 = 
for central collisions (where Npar denotes the number of 
participant for a given centrality). In Fig. |4l we also 
show our predictions for the charged-hadron multiplic- 
ity in pp collisions. We should emphasize again that the 
main difference in our formulation for the case of pp and 
AA collisions is only due to the employed different satura- 
tion scale for a proton and a nuclear target Eqs. (I7I6P (see 
Fig.[T]) and the rest of the formulation is the same. This is 
in accordance with the notion of universality of the satu- 
ration physics which can be further tested at the LHC. In 
Fig.m we show that the recent data from CMS at y/s = 7 
TeV remarkably confirms our predictions 3] for pp colli- 
sions. Our predictions for dNAA/drj at midrapidity rj = 
for central Pb-Pb collisions (for S < 3.7 fm or approxi- 
mately — 6% centrality bin) for ^/s = 2.75 and 5.5 TeV 
are 1152 ±81 and 1314 ± 92, respectively. Similar results 
was also suggested in Ref. [2J| based on the numerical so- 
lution of the BK equation with running strong-coupling 
with the model assumption that the dipole amplitude 
does not depend on the impact-parameter. However, in 



this paper [2J|j the relation between the unintegrated 
gluon-density and the forward dipole-nucleus amphtude 
in the fcr-factorization was taken a simple Fourier trans- 
formation instead of Eqs. ()2|3|) . The predictions of other 
approaches at the LHC can be found in Ref. [25]. 

In Fig. [5l we show [2 / Npar)dN aa / drj as a function of 
number of participant for RHIC and the LHC energies at 
midrapidity. 

Notice that the KLN approach similar to here is also 
based on the /cT-factorization. Both approaches describe 
the RHIC data but provide rather different predictions 
at the LHC. In Fig. HI we also show the KLN predic- 
tions for both pp and AA collisions. Obviously, the KLN 
predictions underestimated the pp multiplicity data at 
the LHC \/s = 7 TeV (in contrast to our predictions) 
while it overestimates the multiplicity for AA collisions 
compared to our predictions. The main differences be- 
tween our approach and the KLN one is that we used 
explicitly the impact-parameter dependent form of the 
/cT-fa-ctorization and employed the correct relation be- 
tween the unintegrated gluon density and the forward 
dipole-nucleon amplitude Eqs. ()2I3|) . Then we employed 
an impact-parameter dependent saturation model which 
gives a good description of HERA data at small- a;. In this 
sense, we did not have any freedom to model the satura- 
tion dipole-proton amplitude at RHIC or the LHC. The 
key difference between our approach and the KLN one 
is the fact that we rely on Eq. (jS]) to determine the sat- 
uration scale for nucleus while in the KLN approach [9] 
the energy dependence was used to obtain the saturation 
scale for nucleus at the LHC from the one at RHIC. 



IV. CONCLUSIONS 



We developed a saturation approach, based on the 
CGC theory, which describes DIS data at HERA [H], 
hadron production at the LHC [3| in pp collisions and at 
RHIC in A A collisions (this paper). In this approach we 
treat in the same way proton and nuclear target consid- 
ering that the only difference between them is the value 
of the saturation scale. The nonperturbative stage of the 
jet-hadronization was described in the same scheme for 
both proton and nuclear target, minimizing the uncer- 
tainties in prediction that could stem from this stage. We 
provided here various quantitative predictions including 
the rapidity, centrality and energy dependencies of the 
inclusive charged-hadron multiplicity in AA collisions at 
the LHC. We believe that our predictions for nucleus- 
nucleus collisions at the LHC will be a crucial test of the 
CGC approach. 
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